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Abstract 

We demonstrate that a tensor product structure could be obtained by introducing pseudorandom 
phase sequences into classical fields with two orthogonal modes. Using classical fields modulated 
with pseudorandom phase sequences, we discuss efficient simulation of several typical quantum 
states, including product state, Bell states, GHZ state, and W state. By performing quadrature 
demodulation scheme, we could obtain the mode status matrix of the simulating classical fields, 
based on which we propose a sequence permutation mechanism to reconstruct the simulated quan- 
tum states. The research on classical simulation of quantum states is important, for it not only 
enables potential practical applications in quantum computation, but also provides useful insights 
into fundamental concepts of quantum mechanics. 
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I. INTRODUCTION 



The classical simulation of quantum systems, especially of quantum entanglement has 
been under investigation for a long time [TH3]. In addition to potential practical applications 
in quantum computation, research on classical simulation systems could help understand 
some fundamental concepts in quantum mechanics. However, it has been pointed out by 
several authors that classical simulation of quantum systems exhibit exponentially scaling 
of physical resources with the number of quantum particles [3j H]. In Ref. [3], an optical 
analogy of quantum systems is introduced, in which the number of light beams and optical 
components required grows exponentially with the number of cebits. In Ref. [5], a classical 
protocol to efficiently simulate any pure-state quantum computation is presented, yet the 
amount of entanglement involved is restricted. In Ref. [1], it is elucidated that in classical 
theory, the state space of a composite system is the Cartesian product of subsystems, whereas 
in quantum theory it is the tensor product. This essential distinction between Cartesian and 
tensor products is precisely the phenomenon of quantum entanglement, and viewed as the 
origin of the limitation of classical simulation of quantum systems. 

In wireless and optical communications, orthogonal pseudorandom sequences have been 
widely applied to Code Division Multiple Access (CDMA) communication technology as a 
way to distinguish different users [7J . A set of pseudorandom sequences is generated by 
using a shift register guided by a Galois field GF(p), that satisfies orthogonal, closure and 
balance properties [Tj. In Phase Shift Keying (PSK) communication systems, pseudoran- 
dom sequences are used to modulate the phase of the electromagnetic/optical wave, where 
a pseudorandom sequence is mapped to a pseudorandom phase sequence (PPS) values in 
{0, 2n/p, ■ ■ ■ ,2k (p — 1) /p}. Guaranteed by the orthogonal property of the PPS, different 
electromagnetic/optical waves could transmit in one communication channel simultaneously 
with no crosstalk, and could be easily distinguished by implementing a quadrature demod- 
ulation measurement [5]. 

In this paper, by introducing the PPSs into classical fields, we explore an efficient simula- 
tion of quantum states based on classical fields with two orthogonal modes. We demonstrate 
that n classical fields modulated with n different PPSs can constitute a 2™-dimensional 
Hilbert space that contains tensor product structure, which is similar with quantum sys- 
tems. In Ref. [H], the efficient classical simulation of Bell states and GHZ state has been 
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introduced and both correlation analysis and von Neumann entropy have been applied to 
characterize the simulation. In this paper, by performing quadrature demodulation scheme, 
we could obtain the mode status matrix of the simulating classical fields, based on which we 
propose a sequence permutation mechanism to reconstruct the simulated quantum states. 
Besides, classical simulation of some other typical quantum states is discussed, including 
product state and W state. We generalize our simulation and discuss the efficiency of our 
simulation in the final end. 

The paper is organized as follows: In Section [TTj we introduce some preparing knowledge 



needed later in this paper. In Section III , the existence of the tensor product structure in our 
simulation is demonstrated and the classical simulation of several typical quantum states is 



analyzed. In Section IV, a generalization of our simulation is proposed and the efficiency of 
the simulation is discussed. Finally, we summarize our conclusions in Section [V] 



II. PREPARING KNOWLEDGE 



In this section, we introduce some notation and basic results needed later in this paper. 
We first introduce pseudorandom sequences and their properties. Then we discuss the 
similarities between classical field and single-particle quantum states. Finally, we introduce 
the scheme of modulation and demodulation on classical fields with PPSs. 



A. Pseudorandom sequences and their properties 

As far as we know, orthogonal pseudorandom sequences have been widely applied to 
CDMA communication technology as a way to distinguish different users [6j [7j. A set of 
pseudorandom sequences is generated by using a shift register guided by a Galois Field 
GF(p), that satisfies orthogonal, closure and balance properties. The orthogonal property 
ensures that sequences of the set are independent and distinguished each other with an 
excellent correlation property. The closure property ensures that any linear combination of 
the sequences remains in the same set. The balance property ensures that the occurrence 
rate of all non-zero-element equals with each other, and the the number of zero-elements is 
exactly one less than the other elements. 

One famous generator of pseudorandom sequences is Linear Feedback Shift Register 
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(LFSR), which could produce a maximal period sequence, called m-sequence [TJ. We con- 
sider an m-sequence of period N — 1 (N = p s ) generated by a primitive polynomial of degree 
s over GF(p). Since the correlation between different shifts of an m-sequence is almost zero, 
they can be used as different codes with their excellent correlation property. In this regard, 
the set of N — 1 m-sequences of length N — 1 could be obtained by cyclic shifting of a single 
m-sequence. 

In this paper, we employ pseudorandom phase sequences (PPSs) with 4-ary phase 
shift modulation, which is a well-known modulation format in wireless and optical com- 
munications, including Orthogonal Quadrature Phase Shift Keying (O-QPSK) and Mini- 
mum Phase Shift Keying (MSK) [9]. We first propose a scheme to generate a PPS set 

over GF(4) [TU] . A^ is an all-0 sequence and other sequences 
can be generated by using the method as follows: 

(1) given a primitive polynomial of degree s over GF(4), a base sequence of a length 4 s — 1 
is generated by using LFSR; 

(2) other sequences are obtained by cyclic shifting of the base sequence; 

(3) by adding a zero-element to the end of each sequence, the occurrence rates of all 
elements in all sequences are equal with each other; 

(4) mapping the elements of the sequences to [0,27r]: mapping 0, 1 mapping tt/2, 2 
mapping n, and 3 mapping 37r/2. 

Further, we define a map / : A — > e lX on the set of S, and obtain a new sequence set 
Q = (pU) = e* A(j) , j — 0, . . . , N — 1 1. According to the properties of m-sequence, we 

can obtain following properties of the set Q, (1) the closure property: the product of any 
sequences remains in the same set; (2) the balance property: in exception to any 
sequence of the set fl satisfy 

N N t c\ \ 

k=l k=l 

(3) the orthogonal property: any two sequences satisfy the following normalized correlation 

1 N 

E(<P®,<P®) = ^E^VP (2) 




In fact, the map / corresponds to the modulation of PPSs of Q on classical fields. Ac- 
cording to the properties above, the classical fields modulated with different PPSs become 
independent and distinguishable. 



B. Similarities between classical field and single-particle quantum states 

We note the similarities between Maxwell equation and Schrodinger equation. In fact, 
some properties of quantum information are wave properties, where the wave need not be a 
quantum wave [3]. Analogous to quantum states, classical fields also obey a superposition 
principle, and could be transformed to any superposition state by unitary transformations. 
Those analogous properties made the simulation of quantum states using polarization or 
transverse modes of classical fields possible [TTHT3] . 

We first consider two orthogonal modes (polarization or transverse), which are denoted 
by |0) and |1) respectively, as the classical simulation of quantum bits (qubits) |0) and |1) 

mm- 

io)= i ; i i) =r i (3) 



Thus, any quantum state \<p) = a\0) + /3\1) can be simulated by a corresponding classical 
mode superposition field, as follows 

|V0 =a\0)+p\l), \a\ 2 + |/3| 2 = 1, (a, $ G C) . (4) 

Obviously, all the mode superposition fields could span a Hilbert space, where we can perform 
unitary transformations to transform the mode state. For example, the unitary transforma- 
tion U (x, 0) is defined 

• I a, ■ a\ / COSY — ie l6 SmY \ , 

U(x,9) = e lx{a * cos0+a y sme) = , (5) 

\ ie~ %e sin x cos X I 

where a x , a y are Pauli matrices. The modes |0) and |1) can be transformed to mode 
superposition by using U (x$) ; respectively, as follows 

U(x,e)\0) = cos X |0)+^sin X |l), (6) 
U(x,9)\l) = cos X |l)-^sin X |0). 
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FIG. 1: The PPS encoding scheme for one input field, where PNG denotes the pseudorandom 
number generator and PM denotes the phase modulator. 



Now, we consider some devices with one input and two outputs, such as beam or mode 
splitters, which split one input field \ipi n ) — a |0) + f5 |1) into two output fields 



and 



ip^t ) • F° r the case of beam splitters, the output fields are 



and 



Yout 



- C a («|0) + /3e^|l)) 

Cb (a |0) + f5e l ^ h |1)) with an arbitrary power ratio |C a | 2 : | C5 1 2 between the 
output beams, where <fi a ,b are the additional phases due to the splitter. For the case of mode 



Yout 



splitters, the output fields are 



Yout 



Yout 



p e i4>b w here 4> a ,b are also 



= ae^ a |0) and 

the additional phases. Conversely, the devices can act as beam or mode combiners in which 
beams or modes from two inputs are combined into one output. 



C. Modulation and demodulation on classical fields with pseudorandom phase 
sequences 

We first consider the modulation process on a classical field with a PPS. Similar to O- 
QPSK system, chosen a PPS A™ in the set of H, the phase of the field could be modulated 
by a phase modulator (PM) that controlled by a pseudorandom number generator (PNG), 
the scheme is shown in Fig. 1. If the input is a single-mode field, it could be transformed 
to mode superposition by performing a unitary transformation after the modulation. 

Then we consider the quadrature demodulation process of a modulated classical field. 
Quadrature demodulation is a coherent detection process that allows the simultaneous mea- 
surement of conjugate quadrature components via homodyning the emerging beams with 
the input and reference fields by using a balanced beam splitting, where the reference field 
is modulated with a PPS \^ r \ The differenced signals of two output detectors are then 
integrated and sampled to yield the decision variable. We can express the demodulation 
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FIG. 2: The PPS quadrature demodulation scheme for one input field, 
process in mathematical form 

i N 

/(AW, A W) = ^£cos(a«-AM) ( 7 ) 

fc=i 




where Ajj/,AjjT are the PPSs of the input and reference fields respectively. The output 
decision variable is 1 if and only if A^ , are equal; otherwise the output decision variable 
is 0. The results are guaranteed by the properties of PPSs. If the input is a single-mode field, 
the scheme shown in Fig. 2 is employed to perform quadrature demodulation. Otherwise 
the scheme shown in Fig. 3 is used, in which the input field = e lX<l) (a* |0) + |1)) is 
first splitted into two fields «je iAW |0) and f3 i e' lXM |1), and two coherent detection processes 
are then performed on the two fields respectively. Noteworthily, the modes of the reference 
fields must be consistent with the two output fields. Thus there are two output decision 
variables a and (3, which correspond to the modes |0) and |1), respectively. We define y5c, pj 
as the mode status. Besides the quadrature demodulation, we can also easily measure the 
amplitudes of the fields |a,| , after mode spitting in the scheme. 



III. CLASSICAL SIMULATION OF MULTIPARTICLE QUANTUM STATES 

We discuss simulation of multiparticle quantum states using classical fields modulated 
with PPSs in this section. We first demonstrate that n classical fields modulated with 
n different PPSs could constitute a similar 2 n -dimensional Hilbert space that contains a 
tensor product structure. Besides, by performing quadrature demodulation scheme, we 
could obtain the mode status matrix of the simulating classical fields, based on which we 
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FIG. 3: The PPS quadrature demodulation scheme for one field with two orthogonal modes, where 
the gray block denote the mode splitter. 

propose a sequence permutation mechanism to reconstruct the simulated quantum states. 
The classical simulation of some typical quantum states is discussed, including product state, 
Bell states, GHZ state and W state. 



A. Classical fields modulated with pseudorandom phase sequences and their ten- 
sor product structure 

For convenience, here we consider two classical fields modulated with PPSs and their 
tensor product structure. Chosen two PPSs of A^ and A^ from the set H, any two fields 
modulated with the PPSs can be expressed as follows, 

m = e iA(a) (a a |0) + /3 a |l)), (8) 
\i? h ) = e iA(6) (a b |0) + /3 6 |l)). 



We define the inner product of two fields \ip a ) and |V>&) as follows, 

N 

X 



N 

= ^J2 e'^'^ + PM ■ (9) 



k=l 



According to the properties of the PPSs, we can easily obtain 



U'„h,) = { 1,a b , (io) 

0, a b 
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which shows that two fields modulated with two different PPSs are orthogonal. The orthog- 
onal property supports to construct the tensor product structure of multiple fields. 

Assume two Hilbert spaces w and v spanned by the states \if) a ) and \ipb), any linear 
combinations of the elements in the direct product space of w <g) v remain in the same 
space. We define the two orthogonal modes modulated with the PPS as the orthonor- 
mal bases of the space of w } expressed as |0 a ) = e lA(a> |0) and |l a ) = e tX(a) |1). Using 
the same notion, the orthonormal bases of the space of v are expressed as |0ft) = e lX(b) |0) 
and | If,) = e iA<6> |1). The four orthonormal bases are thus independent and distinguish- 
able. Then the orthonormal bases of the direct product space of w £g> v can be expressed as 
{|0 a ) <S> | Oft) , |0 a ) <8) |lft) , |l a ) <g) | Oft) , |l a ) <g) | lft)}. Further, we can obtain the following tensor 
product properties, 

(1) for any scalar z, the elements \ip a ) , \ipb) in the spaces of w and v, respectively, satisfy 

Z ® 1^)) = (z ® \A) = |Va) ® (z ; (ii) 

(2) in the space of w (g) v, the direct product of the combinations of elements is equal to 
the combination of the direct products of elements, 

(MO + Wa)) ® (1^) + 1^)) = M ® W + ® 1^) (12) 

+ (s IV%) + ® IV£) • 

Using the same notion, we can construct a 2 ra -dimensional direct product space of |?/>)®" = 
IV'i)®- • •^iV'n) by using the mode superposition of n classical fields l^i) , • • • , \ip n ) modulated 
with n PPSs. 

Quantum entanglement is only defined for Hilbert spaces that have a rigorous tensor 
product structure in terms of subsystems. Here we construct a similar structure of multiple 
classical fields, which is the basis of efficient classical simulation of quantum entanglement. 



B. Reconstruction of quantum states based on the simulating classical fields 



We have discussed quadrature demodulation process in Sec. |II C| Here we discuss how to 
reconstruct quantum state based on the simulating classical fields with the help of quadrature 
demodulation. 
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First, we consider the general form of n classical fields modulated with PPSs 
(A^ 1 -*, . . . , A*™)} chosen from the set H, 



|Vi) = £T (e lAW + ... + e lAW )|0)+(e^ + ... + e lAW )|l) , (13) 



= ^ [(e aW + + e iA(/) ) |0) + + ... + e" w j |1) 

where Cj are the normalized coefficients and a, . . . ,h are the sequence numbers. It is note- 
worthy that although multiple PPSs are superimposed on both modes of the fields, all of the 
PPSs could be demodulated and discriminated by performing the quadrature demodulation 



introduced in Sec. II C . which has already been verified by many actual communication 
systems [9]. 

Now we propose a scheme, as shown in Fig. 4, to perform the quadrature demodulation 



introduced in Sec. II C In the scheme, quadrature demodulations are performed on each 
field, in which the reference PPSs are ergodic on {A™, . . . , A^"H. Thus a mode status 
matrix M (ai,pf\, as shown in Fig. 5, could be obtained by performing n quadrature 
demodulations on the n classical fields. Of the matrix M (a^, pfj , each element is the mode 
status of the ith classical field when the reference PPS is A^. The element takes one of four 
possible values: (1, 0) , (0, 1) , (1, 1) or 0, denote that the PPS A^ is modulated on mode |0) of 
the ith classical field, on mode |1), on both |0) and |1), on neither |0) nor |1), respectively. It 
is noteworthy that different modulation of the n classical fields correspond to different mode 
status matrixes, and vice versa. Thus we obtain a one-to-one correspondence relationship 
between the n classical fields and the mode status matrix. Besides, further discussion will 
show that structure of quantum states and quantum entanglement could be revealed in the 
mode status matrix, which means that a correspondence could also be obtained between the 
mode status matrix and quantum states. Thus we treat the mode status matrix as a bridge 
to connect the simulating fields and the quantum states. 

Now we focus on the correspondence between the mode status matrix and quantum states, 
and consider how to reconstruct quantum states based on the mode status matrix. We first 
transform the matrix M ( ad, m ) to a block diagonal matrix by permuting the fields and the 



sequences, namely the rows and columns in the matrix, respectively, and obtain a matrix 
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expressed as similar to 



/ Ml 



\ 



Ml Ml 
M\ Ml 



M 



M\ Ml Ml 
Ml Ml Ml 
Ml Ml Ml 



(14) 



V 



The structure of the matrix could clearly reveal the structure of the simulated quantum 
state. Each diagonal block of the matrix denotes one unreduced subsystem of the simulated 
state |\&). Thus the state \$>) could be expressed as the direct product of the unreduced 
states |^ 6 ), as follows 



where m denotes the number of the matrix blocks. 

Now we consider how to reconstruct each \^ b ) based on each submatrix, respectively. As 
each submatrix is irreducible, each \ty b ) corresponds to an entangled state. More important, 
different entanglement structures correspond to different structures of the unreduced subma- 
trix. Thus a correspondence relationship could be obtained between quantum entanglement 
and the unreduced mode status matrix. In order to reconstruct the quantum entanglement 
state, an ergodic ensemble of PPSs is required to obtain all possible base states. Thus we 
propose a sequence permutation mechanism to reconstruct each \^ b ) based on each subma- 
trix, which is one of the simplest mechanisms for sequence ergodic ensemble. Assumed \^b) 
contains / classical fields with / PPSs, namely the corresponding submatrix contains / rows 
and I columns, the sequence permutation is arranged as 



R, = {AW, A< 2 ), . . . , A«} , R 2 = {A( 2 \ A (3) , . . . , A«} , . . . , R t = {A« A«, . . . , A ( ' _1) } . (16) 



In the mechanism, each R r corresponds to one selection from the unreduced submatrix. 
We obtain a direct product of I items for each R r , and the simulated quantum state is the 
superposition of the / product items. Therefore we obtain 



rn 



*)=ni*»)> 



(15) 



6=1 




(17) 
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FIG. 4: The PPS quadrature demodulation scheme for multiple input fields, where the M block is 
shown in Fig. 3. 

a i r i P% r ) i s the m ode status obtained from the 

submatrix M (&{, pjj , where j [R l r ] denotes the sequence number of the zth sequence in R r . 

It is noteworthy that the mechanism we proposed above is one of the feasible ways 
to reconstruct the simulated quantum state based on the unreduced mode status matrix. 
Other mechanisms might also work, as long as a sequence ergodic ensemble is obtained in 
the mechanism. The sequence permutation mechanism above could successfully reconstruct 
many quantum states, including the product states, Bell states, GHZ states and W states. 
We will discuss the related contents in next subsection. 



C. Classical simulation of several typical quantum states 

In this subsection, we discuss classical simulation of several typical quantum states, in- 
cluding product state, Bell states, GHZ state and W state. 
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FIG. 5: The mode status matrix related to the fields and PPSs, where M? is the element of 
M (del, j3j) for the ith classical field and the reference PPS 



1. Classical simulation of product state 



First, we discuss classical simulation of n quantum product state. The simulation fields 
are shown as follows 



5 a« 



(|o) + |i)) 



(18) 



v/2 



(|0) + |1)) 



By employing the scheme as shown in Fig. 4, we obtain the mode status matrix 

/ (i,i) o N 



\ 



(19) 



o (M) y 

which demonstrates that each classical field is the superposition of two orthogonal modes 



and no entanglement is involved. According to Eq. (15), we obtain 

1 



I*) 



2 n/2 
1 



(|0) + |1))®...<8>(|0) + |1)) 
(|0...0) + |0...1) + ... + |i...i)) 



(20) 



where \q 1 ...q n ) = 



2 n/2 

S> \q n ) , (ft = or 1) 
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2. Classical simulation of Bell states 



Now we discuss classical simulation of one of the four Bell states = 
^ (|0 a ) |0 b ) + |l a ) | U)), which contains two classical fields as follows 

M = ^(e^ ) |0) + e^ ) |l)), (21) 
M = ^(e iA(b) |0) + e- W |l)). 
By employing the scheme as shown in Fig. 4, we obtain the mode status matrix 

M 1*1,01) = I [1,0) ^'^ I (22) 



(0,1) (1,0) 

We note that in this case, the mode status matrix is irreducible, which corresponds to an 
entanglement state. According to the sequence permutation mechanism, we obtain that 
Ri = {A^,A^} and R 2 = {\( b \\^}. Based on the mode status matrix, for the selection 
of Ri, we obtain |0) <8> |0); for the selection of R 2 , we obtain |1) ® |1). If we randomly choose 
one selection between R x and R 2 , we could randomly obtain one result between |0) ® |0) and 
|1) <g) |1), which is similar with the case of quantum measurement for the Bell state \^ + ). 
We could reconstruct the state based on the mode status matrix 

2 



* + ) = j= E i°) + ft™ i 1 )) ® io) + ii) 



^(|00) + |11)) 



(23) 



In quantum mechanics, another Bell state |$ + ) could be obtained from by per- 

forming the unitary transformation a x : |0) -H- |1) on one of the particles. Using the same 
method, we perform an unitary transformation on \ipb) to flip its modes |0) -H- |1). Thus we 
obtain two classical fields as follows 

m = j=(e^\0) + e^\l)), (24) 



V2 

\A) = j=(e* w \l) + e» w \oj). 



V2 

By employing the scheme as shown in Fig. 4, we obtain the mode status matrix 

M^H^^I. (25) 



(1,0) (0,1) 



14 



According to the sequence permutation mechanism, here we obtain Ri = {A^, A^} and 
R% = {A^, A^} again. As the mode status matrix is different, for Ri, the result turns to 
be |0) ® |1); for R2, we obtain |1) (g) |0). If we randomly choose one selection between R\ and 
i?2, we could also randomly obtain one result between |0) Cg> |1) and |1) <S> |0), which is similar 
with the case of quantum measurement for the Bell state |$ + ). We could reconstruct the 
state 

2 r 



|$+) 



1 

71 



V2 



E 

r=l 



(|01) + |10)). 



ajwi |o) + /3^]|i))®(« b Kj |o) + /3; 



ir 



(26) 



For other two Bell states \^~) and |$~), they could be obtained from \^ + ) and |$ + ) by 
using a tt phase transformation. However, they could not be distinguished from \^ + ) and 
|$ + ) unless an orthogonal projection measurement is performed [8]. 



3. Classical simulation of GHZ state 



For tripartite systems there are only two different classes of genuine tripartite entangle- 
ment, the GHZ class and the W class [HI [T3j. First we discuss the classical simulation of 
GHZ state \^ghz) — 771 (|0 a ) |0fc) |0 C ) + |l a ) 1 1^} |l c )), which contains three classical fields as 
follows 

1 



IV*) 



V2 



|0) + e a(6) |l 



(27) 



1 V b) |0) + e^ (c) |l 



V2 
1 

71 



e ^ c) \0) + e aw \l) 



Performing the scheme as shown in Fig. 4, we obtain the mode status matrix 

( (1,0) (0,1) ^ 

(4$ 



M 



(28) 



(1,0) (0,1) 
^(0,1) (1,0) ) 

According to the sequence permutation mechanism, we obtain that R\ = {X^ a \ X^ b \ A^}, 
R 2 = {A (6) , A (c) , A (a) } and R 3 = {A (c) , A (a) , A (fe) }. Based on the mode status matrix, for the 
selection of Ri, we obtain |0) ® |0) ® |0); for the selection of R 2 , we obtain |1) <g> |1) <g> |1); 
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for the selection of i? 3 , we obtain nothing. Thus we could reconstruct the state based on 
the mode status matrix 



V2^ 

v r=l 
1 

V2 



\^hz) = ± £ fe™ 1°) + & m I 1 )) ® 1°) + S K1 1 1 )) ® (<^ ] 10) + |i) 

(|000) + |111)). (29) 



4- Classical simulation of W state 

Then we discuss the classical simulation of W state, 

\V W ) = -i= (|1 ) |0 b ) |0 C ) + |0 o ) |1 6 ) |0 C ) + |0 a ) |0 b ) |l c )) , 



(30) 



which contains three classical fields as follows 



\A) = ^(e a(a) |l) + e^ 



|0) + e tA(c) |0)), 



(31) 



1 



[e— |l) + e a(6) |0) + e lAW |0)), 



|l) + e iA(b) |0) + e a(c) |0)), 



It is noteworthy that the three classical fields could be produced from one single field by 
using two beam splitters, which is quite similar with the generation of W state in quantum 
mechanics. Performing the same scheme, we obtain the mode status matrix 



( (0,1) (1,0) (1,0)^ 



(32) 



(0,1) (1,0) (1,0) 
^(0,1) (1,0) (1,0) J 

According to the sequence permutation mechanism, we obtain that R\ = {A^,A^,A^}, 
R 2 = {A^, A (c) , A (a) } and R 3 = {\ {c \ A (a) , A (b) } again. Based on the mode status matrix, we 
obtain |l)<g> |0) <g> |0), |0)<g> |0)<g> |1), |0) <g) |1) ® |0) for the selection of R ± , R 2 , R 3 , respectively. 
We find an interesting fact that when we obtain the state |1) of the first field, R\ must be 
selected, thus only the |0) ® |0) state could be obtained from the other two fields; otherwise 
when we obtain the state |0) of the first field, the selection could be R 2 or i? 3 , thus the 
state of |0) ® |1) + |1) <8> |0) could be obtained from the other two fields. This fact is quite 
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similar with the case of quantum measurement and the collapse phenomenon for W state in 
quantum mechanics. We could reconstruct the state based on the mode status matrix 

(aj*a |o) + M K] |i)) ® (4 R ' ] 10) + K [Rl] ® |0) + |1)) 

= -L(|100) + |010) + |001)). (33) 

IV. GENERALIZATION OF THE SIMULATION AND SOME DISCUSSIONS ON 
EFFICIENCY 

In this section, we first propose a generalization of our simulation to the case of an 
arbitrary quantum state, then the efficiency of our simulation is discussed. Besides, for 
better understanding our simulation, the PPSs in the cases of modulating two and three 
classical fields is illustrated. 



I* 



w 



1 

73 



A. Generalization of the simulation 



Here we propose a generalization of our simulation to the case of an arbitrary quantum 



state. In Sec. Ill C classical simulation of several typical quantum states has been intro- 
duced, including product state, Bell states, GHZ state, and W state, based on which we 
assume that any quantum state of n — 1 particles could be successfully simulated. Here 
we consider the classical simulation of an arbitrary quantum state of n particles, \^ n )- We 
first transform \^/ n ) to an equivalent expression, \ty n ) = |0) n + |0 n -i) |l) n , where 

|$ n _i) and |G n _i) denote two quantum states of n — 1 particles. According to the as- 
sumption, we obtain that the quantum state |$ n -i) could be simulated by n — 1 classical 
fields \(fi) , . . . , \<Pn-i) modulated with n — 1 PPSs {A^, . . . , A^™ -1 )}, and the quantum state 
|0 n _i) could be simulated by n — 1 classical fields \&\) , . . . , |i? n _i) modulated with n — 1 
PPSs {A« ...,A(- 2 ),AW}. 

Further, we consider the classical simulation of the quantum state \^f n ), which is the 
superposition of the classical simulation of |3> n -i) ® |0) and the classical simulation of 

|e„_ x > <s> |i> 



= ^r[ki) + h?i)], (34) 
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[\<P n -l) + |#n-l)] 



l r 




|0) + e 



,iA ( " _1) 



1) 



C n - 



Thus by only adding one classical field and one PPS, we obtain the classical simulation 
of the quantum state |\P n ) based on the assumed classical simulation of |$„_i) and |O n _i). 
Using the principle of induction, we could provide classical simulation of any quantum state. 
Therefore we successfully generalize our simulation. 

B. Efficiency of the classical simulation 

For a long time, researchers have used classical fields to simulate quantum states and 
quantum computation. In these researches, multiple qubits are distinguished by different 
degrees of freedom, such as optical modes or space positions. However, as no tensor product 
structure is obtained in these classical simulations, each quantum base state needs inde- 
pendent degree of freedom to simulate. In quantum mechanics, the number of quantum 
base states grows exponentially with the number of quantum particles. Therefore the classi- 
cal simulation require resources that also grow exponentially with the number of simulated 
quantum particles, which is not efficient. 

In this paper, we utilize the properties of PPSs to distinguish classical fields that are even 
overlapped in same space and time. A 2 n -dimensional Hilbert space which contains tensor 
product structure is spanned by n classical fields modulated with PPSs. In our scheme, the 
resources required are classical fields modulated with PPSs instead of optical/space modes. 
One classical field modulated with one PPS can simulate one quantum particle. It means 
that the amount of classical fields and PPSs grows linearly with the number of quantum 
particles. According to the m-sequence theory, the number of PPSs in the set S equals 
to the length of sequences, which means that the time resource (the length of sequence) 
required also grows linearly with the number of the particles. Based on the analysis above, 
we conclude that one can efficiently simulate quantum entanglement with linearly growing 
resources by using our scheme. 
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C. Illustration of some pseudorandom phase sequences 



For better understanding our scheme, the PPSs in the cases of modulating two and three 
classical fields is illustrated below. Using the method mentioned in section [II], an m-sequence 
of length 4 2 — 1 is generated by a primitive polynomial of the lowest degree over GF(4), which 



is 



120332301131022 



. Further we obtain a group that includes 16 PPSs 
of length 16: {A^ -*, . . . , A^ 15 * 1 }, of which in exception to \(°\ all PPSs are independent and 
could be used to modulate classical fields to simulate quantum states of up to 15 particles. 
We could choose any two PPSs from the group for the simulation of two particles quantum 
state, for example, 



A (.) 



1203323011310220 
2120332301131020 



xtt/2, 
x tt/2. 



(35) 



And we could also choose any three PPSs from the group for simulation of three particles 
quantum state, for example, 



A («) 

A <6) 
A (c) 



1203323011310220 
2120332301131020 
2212033230113100 



x tt/2, 
xtt/2, 
x tt/2. 



(36) 



Using classical fields modulated the PPSs given above, we can efficiently simulate any quan- 
tum state of two or three particles. 



V. CONCLUSIONS 



In this paper, we have discussed a new scheme to simulate quantum states by using 
classical fields modulated with pseudorandom phase sequences. We first demonstrated that 
n classical fields modulated with n different PPSs can constitute a 2 n -dimensional Hilbert 
space that contains tensor product structure, which is similar with quantum systems. Fur- 
ther, by performing quadrature demodulation scheme, we obtained the mode status matrix 
of the simulating classical fields, based on which we proposed a sequence permutation mech- 
anism to reconstruct the simulated quantum states. Besides, classical simulation of several 
typical quantum states was discussed, including product state, Bell states, GHZ state and 
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W state. We generalized our simulation and discussed the efficiency of our simulation fi- 
nally. We conclude that quantum states can be efficiently simulated by using classical fields 
modulated with pseudorandom phase sequences. The research on simulation of quantum 
states may be important, for it not only provides useful insights into fundamental features of 
quantum mechanics, but also yields new insights into quantum computation and quantum 
communication. 
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Fig. 1 The PPS encoding scheme for one input field, where PNG denotes the pseudorandom 
number generator and PM denotes the phase modulator. 

Fig. 2 The PPS quadrature demodulation scheme for one input field. 

Fig. 3 The PPS quadrature demodulation scheme for one field with two orthogonal modes, 
where the gray block denote the mode splitter. 

Fig. 4 The PPS quadrature demodulation scheme for multiple input fields, where the M 
block is shown in Fig. 3. 

Fig. 5 The mode status matrix related to the fields and PPSs, where M\ is the element of 
M (aj, f3f ) for the iih classical field and the reference PPS 
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